On the Levy density function 
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Abstract 

In this paper, we introduce the Levy density function as the hmit of a general- 
ized Mittag-Leffler density function. The fractional integral equation for the generalized 
Mittag-Leffler density function is also given. And the role of the Levy structure in the 
I fractional calculus is described. Finally, a transformation is defined. 
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1 Introduction 

In the following sections, we are going to look into some details related to Levy distribution 
and its role in the fractional calculus. We state some preliminary results. Pochammer symbol 
is defined as 

{b)k = b{b+l)---{b + k-l), (6)o = l, MO. (LI) 



' Definition 1.1 The Euler gamma function T{z) is defined as follows 

o 

(N 



POO 

T{z) = e~^e-P*dt, > 0, 9^(z) > (1.2) 

JO 

g; = lim — z/0,1,2,3,--- . (1.3) 

' n^oo z[Z + 1) ■ ■ ■ [Z + n) 

The representation of the Pochammer symbol in terms of gamma functions is 

(»).^^^ (1.4, 



whenever the gammas exist. 

Lemma 1.1 / Stirling asymptotic formula J f^l^ 
For 1^1 — )■ oo and a a bounded quantity, 



T{Z + a) « (2^)l/2^^+a-l/2g-^^ (-^^5) 



Lemma 1.2 For |7[ ^ oo and k a bounded quantity, 



(7)fc r(7 + k) 

7— >oo ■7''' 7— >oo 1 (7)7 
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/f- function representations and their convergent regions are important. We are not going 
to explain them and the definition of Mittag-Leffler function, its generalized ones and their 
properties are not added here. So for those, [3l [9l [10] have detailed descriptions. 

In the paper [S] , the author used some properties of Mellin transformation and statistical 
methods to investigate the Mittag-Leffler statistical distribution and its generalized ones. In 
[6], he showed the possibility of getting the explicit form of Levy density function. He also 
computed E[vaY~^ from a gamma density ([6] for details). This computation will lead us to 
define a transformation in section 4. In [8J, the following results are stated. 



Lemma 1.3 For > 0,^(7) > 0, 



i/3i^oo r(7) 



rl,l 
1,2 



i^^^)%^uU.)\=i^-A-' (1-7) 



Theorem 1.1 For 3fJ(/3) > 0, 3^(7) > 0, x > 0, a > 0, g > 1, c> 0, 



lim ^, „ , x ^ E( r.. 



a{q-l){xpfy =cx^ l + a(g-l)x' (1. 



Theorem 1 1.11 is related to the pathway model, for x > from where one has Tsallis statistics, 
superstatistics, power law and many others. 

The concept of pathway model is important in the sense that it is a pathway between two 
totally different looking systems. It can be considered that the behavior of the family of den- 
sity functions tied up by the pathway parameter is similar to that of the collection of germs 
connecting two different points via a path in the complex plane. For details, see [U EJ El E] . 

In this paper, we construct the Levy density function through a pathway parameter. Fi- 
nally we define a transformation of connecting ordinary space and a-fractional (or a-level) 
space. 

2 The Levy density function 

Consider a gamma random variable x with the density function 

9{x) 



^7^7 Ig 'yx 



0, otherwise. 
This density function has some interesting properties. 
1. Laplace transform : L,{s) = (l + 



jy—. for < X < 00, 7 > 

r(7) ' ' (2.1) 



2. Mellin transform : Mg{s) = ^^^^t£_li^i-^ 
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3. Mellin-Barnes integral representation : 



1 / r(7 + s-i) , 



27ri Jl, r(7) 
We consider 

/(^) = ^""-S^i?J.,,,)(-7x"), (2.2) 

which can be obtained by using the technique developed in section 4 with the density function 
in (|2.ip . This is a generalized Mittag-Leffler density function and the function which leads to 
the explicit form of Levy density function. 



1. Laplace transform : 



2. Mellin transform 



Lf{s) = [l + -] . (2.3) 



^/(^) ar(7)r(i - .) ^ " " ^^-^^ 



where < i2e(s) < a < 1. 
3. Mellin-Barnes integral representation : 

1 / r(s)r(7-s) 



27ri r(7)r(a7 — as 

1 f r(7 + ^ - i)r(-^ + i 



^a7-as-l^-s+7^g (2.5) 

7-f+^a;-^cis (2.6) 



27ri Jl> ar(7)r(l - s) 
where < Re{s) < a < 1. 

Firstly, we look into its Laplace transform. When 7 tends to 00 in (j2.3p . (|2.3p goes to 
e~*" which is the Laplace transform of Levy distribution. This means that (|2.2p will lead 
us to the Levy density function. For this, we make use of (|2.6p and lemma 11.21 From the 
representation (12. 6p . we have 

1 I r(7+i^i)r(-^+i) . 

7 at ax as. 



27ri f]^' Q!r(7)r(l — s 
As 7 tends to 00, we get 

1 f r(-^ + i 



^ . r X./. ^ x-'ds, < Re(s) <a<l (2.7) 
ZTTZ J II ai (1 — sj 

which is the Mellin-Barnes integral representation of Levy density function [6]. 

Definition 2.1 Let f(x) have H-function representation and be convergent f3l [W\^. Then f(x) 

is said to be a function with Levy structure if its H-function representation has the factor 
r(--+— ) r(— — +-) 

Q,r(i-s") integrand. In short, ^-pii-s) called Levy structure. 
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In [5], the 2- parameter Weibull density function is defined as follows, 

6bx^-^e-'"''' for < x < oo, 5 > 0, b>0 
0, elsewhere. 

Consider the remaining part of (j2.6p after removing the Levy structure, 

1 /• r(7 + ^ - ^) 



(2.8) 



2iTi 



r(7) 



"7 a + QX '^ds. 



(2.9) 



Then by the residue theorem, it becomes y^'^"'^ i^" ^ where a > 0,7 > 0, which looks 
the same as ^^^p^_^^ when a is replaced by 1. This is a density function and a 2-parameter 
generalized gamma density function, we just list 3 cases in the below table. 

generalized gamma density functions from Mittag-Leffler density functions 



generalized gamma density functions 


Mittag-Leffler density functions 






ct-y'' „a7-lp-7a;° 




aS^ a-q—l—Sx" 





When 7 ^ oo in ()2.9p . the Mellin transform becomes 1, which means 

v7 _ „, r(7 + ^ - ^) . 



POO 

lim / x^ 

7-^oo Jo 



-1 «7 _^a^-l^-^x-^^ ^ ^-j^^ 



-7 



r(7) 7-i>oo r(7) 

So as 7 takes very large value, then the parameter s becomes redundant. 



1. 



(2.10) 



3 The fractional integral equation as an extension of the re- 
action rate model 

In [8], they consider fractional integral equations as extensions of the reaction rate model. 
'^^J^^ = —cN{x), c > =^ N{x) — Nq = —c J N{x)dx : reaction rate model which gives 
N{x) — Nq = — c° qD~°'N{x) and from there 

N{x) - Nof{x) = -c" oD-'^Nix) 

where f{x) is a general integrable function on the finite interval [0,6]. 
By applying Laplace transformation. 



(3.1) 



N{s) - NoF{s) = -c" s"iV(s) 
where F{s) is the Laplace transform of f{x) and N{s) is the Laplace transform of N(x). Then 

Nos''F{s) 



N{s) 



NoF{s)^{-l)''c''h- 



ak 



k=0 



4 



= NoF{s) + NoF{s) ^(-l)'=c"'=S-"'= = NoF{s) - NoF{s) ^^-l)k^^k+a^-ak-a 
k=l k=0 

But note that Y^'^^^i^-i f^ak+a^-ak-a jg ^Yie Laplace transform of ^"-^^(a,^) (-c° 
applying the inverse Laplace transformation, 



N{x) 



k„ak ex 



Nofix) + NoY, p. ,^ , 
^1 r(aA;) Jo 

Nofix) - Noc^ Hx - t)"-ii?(,,„)(-c"(x - tr)f{t)dt. 



]. By 

(3.2) 
(3.3) 



It seems that the kernel of the integral in (j3.3p has Levy structure and that if the kernel 
of an integral operator has Levy structure, it goes well with the function possessing Levy 
structure, in other words, it brings more nice forms and if f{x) has the Levy structure, then 
the solution has the Levy structure, too. Here we give some examples. 

Let f{x) be 1 in ()3.ip . then N{x) = NQEa{—c°^x°') from ()3.3p . This can be written as follows. 
functions without Levy structure 







N{x) 


1 




NoE^i-c^x^) 


X 




A^02;^{a,2)(-C"x-) 














xi^-^ 




iVoX^-'i^(a,^)(-c"x") 






iVox^-'i^[j;)(-c"x-) 



functions with Levy structure 



fix) 


— > 


N{x) 






iVox-ii?(,,„)(-c-x-) 


(a,a) ^ ' 




iVox°-ii?^+^^(-c"x") 



Now, we are in a position to show that the density function (12. 2p has its own integral 
equation. By changing to 7, ()3.ip becomes 



N{x) - Nofix) = -7 oD-'^Nix). 
Then from (|3.3p . the solution is accordingly 

N{x) = Nofix) - iVo7 Pix - (-7(x - t)")/(t)dt 



(3.4) 



Nofix) - iVo7 / (x - t)"-^if 



1,2 



7(x-t)^ 



1(0,1) 

l{0,l){l~ci,a) 



/(i)dt 
(x - t)-^-i^dsfit)dt. 
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If we consider this particular equation, then ^'y'^ Ej^ ^^^{—^x"') becomes the solution 

of 

N{x) - iVox°^-S^i?5;i^)(-7x") = -7 oD-'^Nix), (3.5) 
which is the main generalized Mittag-Leffler density function in section 2. 

4 A lifting from ordinary space to Q;-fractional(a-level) space 

In [6], he shows a process to lift a gamma density to a generalized Mittag-Leffler density 
function by using statistical techniques and it is in Example 14.11 

Example 4.1 Let x be a simple exponential random variable with the density function f[x) = 
e~^ . We attach the Levy structure to E{x'^Y~'^ = f^{x'^)^~^e~''^dx. 
Then 

- ^(.,«)(-- ) = ^/^ ^f(r3^^ (4.1) 

by the residue theorem. 

In this section, we will show through a new process how to lift a function in the ordinary 
space to a corresponding function in the a-fractional(or a-level) space. 

Let f{x) be a function, which does not have Levy structure and lives in the ordinary space 
and h{x) be the Levy density function as a kernel. 

Definition 4.1 Define a transform of f{x) 

^<^>(^' - (^) ((f) ) T' ^' f^'^' 

where a fixed in < a < 1, x^O, f2{x) = xf{x) and f is integrable and continuous on the 
interval. 

Example 4.2 Let x be a exponential random variable with the density function 

-l)'x' 



f{x) = e ^, f2{x) = x^ 

1=0 



Then 



J(/)(x) = hm / r-)h t) ^"""^ dt 



7->oo_/q \X J VVt// t 



,^oo J \t) ^ l\ ' ' ^ t 

•^'J 1=0 



lim / ^i-^y v -^-'^ fxnr^'S^V 

7^oo n ^ ^ ' ^ klT{ak + a7) u 



— ' -lYn--^ , ,.7-1 7^ (7)fc(-l)V(^<' rf^ 
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^' ^ ^ fc!r(a/c + a7) Jo 
(/ix 7 so i/iai i/ie integral exists,) 

v-^nn^-^ V ^ 2-^ fc!r(aA; + a7) r(aA; + 07 — a/ — a + 1) 



7— >-oo 



i=0 
00 



fc=0 



= lim ^ 7^ 



7— >-Cx3 ' I! 

i=0 



r(s)r(7 — s)r(— as + ^7 — a/ — 0)7 



-as+a7— 1 



lim 



I r(7)r(a7 — as)r(--as + 07 — — a + 1) 

(-ly 1 /■r(7 + i-i)r(-^ + i)r(i- 



7—^00 



«=0 



al — 0)7 a~^aX 



r(7)ar(l - s)r(2 -s-al-a) 



-ds 



-ds 



{where < R{s) < a < 1) 

^(-ly 1 / T(-f + l)r(l-s-a/-a).T- 
^ n 2m Jl ar{l - s)r(2 -s-al-a) 



1=0 

r( ' 



-ds 



+ - 

a a' 



aT{l 



has no poles in the interval < R{s) < a < 1, 



, . , r(l — s — al — a) , . . , 

so we proceed with — — ; -by putting 1 — s — al — a 



T{2-s-al- a) 
{-ly 1 / r(/ + 1 - asi)r(-asi)x-°^^i+°^'+°^-i 



-asi, 



i 



U 2m J L T{al + a - asi)T{l - asi) 



ds 



= E 



i/ie integrand can have a pole at and only one pole there since < a < 1 

-1)' r(i + 
Ti r(a/ + a) 

x"-'^(a,a)(-x")- 



1=0 

„a-l 



In this example, we have connected two important functions the exponential function in the 
ordinary space and a-exponential function in the a-fractional space. 

CORRESPONDENCE 



ordinary space 



a-level space 
^ 



X" 



g— oXg— 6a; g— (a+6)x 



^"~^^(a,a)(-(« + &K) 



oFi{;a;-t) 



a—l v^oo 



Z^k=0 (a)fcr(afc+a) ' 



a > 1 



iFo(a; ; -t) 



r(aA:+a) 



a—l ■^->oo (a)fe(— -. l ^ -i 



iFi(a;6; -t) 



'fc=0 (b)fcr(afc+Q:) 
{a),,{b)k{-ir^ 



2-Fi(a,6; c; -t) 



E 



00 

k=0 



(c)fcr(afc+a) 



, a>l,6>l,c>l 



5 Conclusion 



It seems that the Levy structure is important in the area of Fractional Calculus in the following 
sense, 

• fractional equations naturally possess the Levy structure such as 

• functions with Levy structure matches fittingly with functions with Levy structure, for 
example the tables shown in section 3 

From section 3 and 4, it is urged that the theory of fractional calculus be developed 
according to the origin of functions where the origin means the level of functions such as 
a-level. 

For this, see [21 [9]. 

Furthermore, if we use the asymptotic behavior of the gamma function , ax") 

can be approximately used in place of X^^^q ' ii^mely, the coefficient of x can 

be generated by using the limiting process [8]. 
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